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The long-time behaviors of the velocity autocorrelation function (VACF) for granular 
flows based on the optimal velocity fluid model are investigated theoretically. It is found the 
long-time tail of VACF disappears in a granular flow under the gravity. 

§1. Introduction 

The long-time tails in the current autocorrelation functions are one of the im- 
portant characteristics in nonequilibrium physics.DS'EJS It is known that the long 
tails are originated from the correlated collisions which cause anomalous transport 
behaviors such as the system size dependence of the transport coefficients.'^''^ The 
analysis of such behaviors is still an active subject in nonequilibrium physics. 

Although the long-time tail with the time t and the spatial dimension 

d even at equilibrium is well established, the interest in the long-tails in granular 
fluids is rapidly growing.l3''3'''^''I® Most of the arguments focus on the velocity 
autocorrelation function (VACF) C{t) which is defined by 



1 ^ 

Jm^ — ^(^i(0)-^,(t)), (M) 



where Vi[t) is the velocity of i— th particle at time t, and V is the volume of the 
system. The average (• • • ) is taken over many different ensembles. 

So far theories and simulations to characterize granular fluids discussed in the 
case of freely-cooling cas^''^ or the flow under the simple shear.H'' From these 
papers, we find that (i) the long tail for VACF and the time-correlation of the shear 
stress exis but the time-correlation for the heat flux decays exponentially in freely 
cooling casCj^I^ and (ii) there is a cross-over from t~'^/'^ to for sheared 

granular gases.'^ On the other hand, the existence of conventional long-time tails 
in VACF of the granular flow in a hopper has been reported experimently,® while 
there is no corresponding theoretical or numerical argument. 

In this paper, thus, we wish to clarify whether there exists the long tail of VACF 
in the granular flow under the gravity. Our theoretical method is based on the 
classical one developed by Ernst et aLl^'^^^ In the next section, we will introduce 
the model we use. In § 3, we will present the details of the analysis. In § 4, we 
will discuss and conclude our results. In Appendix A, we show the details of the 
calculation of C{t). 
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§2. Model 

In this section, we introduce the model what we are analyzing. Although the 
experiment by Orpe and KudrollP^ is designed for the dense flow in which the in- 
fluence of the interstitial air is small, our theoretical tool based on the method by 
Ernst et al^ might not be adequate for dense flows. On the other hand, the drag 
from the side walls or the interstitial fluids plays a crucial role to achieve a steady 
granular flow. Without the drag, particles are accelerated by the gravity and cannot 
reach a steady state. Here, we introduce a fluid model to describe the granular flow 
under the gravity, which is a couple of equations for the mass conservation and the 
momentum conservation: 



dtn + V ■ (nu) = 0, (2-1) 

T 1 

dtU + u-Vu = Vn + uV'^u + -[U(n)ez - u], (2-2) 

mn T 



where n and u are the number density of grains and the coarse-grained velocity 
fleld, respectively. The parameters T, m, v and r repsectively represent the granular 
temperature, the mass of a particle, the effective kinetic viscosity and the relaxation 
time to the optimal velocity U{n)ez with the unit vector along z— axis. For granular 
particles immersed in fluids, U{n) may be represented h:)^ Ui{l - 4>f/{l + 2(j) + 
lA92(p{l — (j))^) where (j) is the volume fraction and Ui = 2a?{pp — pf)g/9rj is the one- 
particle sedimentation velocity with the radius of the particle a, the gravitational 
acceleration g, the density of the particle pp, the density of the fluid pf and the fluid 
viscosity r/. For traffic flows, we often use a phenomenological function for U{n)W^ 
In the model, we assume that the relaxation of the temperature is much faster than 
the other hyrodynamic variables. This fast relaxation has been confirmed at least in 
the case of freely cooling granular gases.l^ This model was original introduced as a 
one-dimensional traffic flow model by Karner and Konhauser.'Sf Similar models have 
been used for the description of the fluidized granular particles under the influence of 
air.13'13'13 Later, we have recognized that traffic fiows and granular flows belong 
to the same universality class.'^''^ It should be noted that the models for fluidized 
bedPJ'1SJ'13 contain the pressure which depends on the density, while our model 
adopts the equation of the state for the ideal gas. Although our model ignores 
the finite density effects in the pressure, the pressure includes the most important 
contribution to keep a system uniform. On the other hand, the density dependence 
of the optimal velocity U (n) is essentially important to describe the emergence of the 
density waves. Therefore, we expect that the essence of hyrodynamic behaviors of 
granular fiows is independent of the choice of the basic equation. 

Let us call the model in eqs. (j2-ip and (j2-2p the optimal velocity fiuid model, 
because the optimal velocity U{n)ez plays important roles. We believe that the 
analysis of the optimal velocity model gives a prototype of granular fiows and traffic 
fiows. 
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§3. Theoretical analysis of VACF 

VACF is characterized by the fluctuations of hydro dynamic variables around the 
uniform flow 

n{r,t)=no, u{r,t) = Uoe^, (3-1) 

where Uq is the abbreviation of U^uq). Ernst et al?- demonstrated that VACF can 
be approximated by 

/r df^ 
dvoUvo)vo ■ J j^Uk{t)P-k{t), (3-2) 

where vq is the relative one-particle velocity to the optimal velocity Uoe^, fo{vo) 
is the one-particle velocity distribution function, Uk.{t) and Ph{t) are respectively 
the Fourier transforms of u[r, t) — UqEz and P(r, t) with the probability density of 
a tracer particle P{r,t). From eq. (l3-2p . we can evaluate VACF once we know the 
behavior of Uk{t) and Pk{t)- 

In general, it is not easy to solve the set of nonlinear equations (j2-ip and (j2-2p . 
When we restrict our interest, however, to the case of linearly stable situation around 
the uniform flow (j3-ip . we can simplify the argument. Such situations can be achieved 
for dilute flows or dense flows. For simplicity, we only consider such linearly stable 
situations. 

The probability density should obey the diffusion equation with the drift velocity: 
{dt + Uodz)P{r,t) = DV'^P{r,t), where D is the diffusion constant. Its solution is 
immediately obtained as 

Pkit) = Pfc(O) exp[-ikM - Dkh]. (3-3) 

On the other hand, it is straightforward to obtain d+1 eigenvalues from the linearzed 
equations of (|2-ip and (j2-2p . which are d — 1 degenerated (T_l = —ik-^UQ — - — vk"^ 



and 

1 

with 



T J \ \T 



2 



^ + 2ikzUo + vkA ±\ (- + 2ikzUo + uk^] - Ax 



(3-4) 



X = k' kiU^ -ikA " " " + Uouk' , (3-5) 

m \ 'T J 

where Uq = dU/d n\yi=yiQ. The most relevant eigenvalue to characterize the slow 
relaxation in the vicinity of /c ^ is cr+. The uniform flow is no longer stable, when 
the real part of cj+ can be positive if the condition Cg > T/m is satisfied, where 

Co = UqUq. 

For simplicity, however, we restrict our interest to the case of weakly stable case 
in which T/m > Cg is always satisfied in this paper. Therefore, the approximate 
solutions of the linearized equations of (|2-ip and (j2-2p are given by 

e°"+* T Tk± 

Uk^{t) ~ [ ik±nk{0) 7 r(fc • Mfc(O)) 

cr+ — (T_ mng m(cj+ — a±) 
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+ \a+ + ik,Uo+ ^ ' ' ° 0/ ' \uk^{Q)l (3-6) 

I cr+ - cr_L J 

rikiO) + -]ik±- itfcx(0 ) 

cr+ — (7_ \ mno t J o'+ — cr_L \ mno r / 

/ k'^T \ 
+ {^^ r + a+ + ik,Uo)ukAO)], (3-7) 

where k± and u^j^ are the wave vector in the horizontal plane and the velocity in 
the horizontal plane, respectively. Here, we stress that this solution is valid when 
exp(— t/r) is negligibly small. If r is large enough, the relaxation term to the optimal 
velocity is negilible and the total energy of the system increases with time. 
Substituting ([331), ([Ml) and ([331) into ([321) we obtain 



where 



Cit) = id-l)C4t) + C,it), (3-8) 
T f dk e^p[{a+ + ikzUo - Dk^)t] 



{2'kY a+ - o-_ 

T{e -kl) + ik,{n^U'jT) 

m(cJ+ — (T_L 



+ a++ikM, (3-9) 



and 



T ( dk exp[(a+ + zfc,[/o-m2)t] [ kjT 

am J (ivrj" cr+ — U- |^ m(c7+ — a±) 

(3-10) 

where we use / dvvfo{v) = 0, / dvv'^fo{v) = nT, Uk.{0) ~ (vo/no), and Pfc(O) ~ 
The evaluation of integrals in eqs. (|3-9p and (j3-10p is straightforward (The details of 
calculation for d = 3 are presented in Appendix). The long time behavior of Cx{t) 
and Cz{t) for d = 3 are respectively given by 

r^T^ e-^o*/(4a) e-^o*/(4a) 

^ "4vr3/2m2(I) + ^)aV2 t5/2 (^-ll) 

and 

96^3/2(^ + Zl)«5/2 i3/2 ^3/2 ' ^"^'^^^ 

where cq = nof/g and a = L' + r(T/m — Cq). From eqs. ()3-lip and ()3-12p . it is obvious 
that Cz{t) has the dominant contribution. Therefore, we obtain 

—c'^t/{Aa) 

C(f):^a(t)oc " ^3^^ , (3-13) 

where the amplitude becomes negative when the system is in the weakly stable region 
of uniform flow. 
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It should be noted that we neglect the fast contribution C± (t) 
arising from a±, which corresponds to the correlation from the shear mode. This 
term becomes relevant when the effects of the gravity and the relaxation to the 
optimal velocity is negligible to recover the conventional tail in the limit of 

r — > oo. In this sense, our analysis presented here is valid for r < a/cg. 

Thus, the long-time behavior of VACF of the granular flow governed by eqs. 
(12-11) and (l2-2p has the exponential damping term as C{t) ~ — e~'^o*/{4a)y'^3/2 
is different from the case without the existence of systematic flows. 

§4. Discussion and Conclusion 

In this paper we have demonstrated the absence of the long-time tail in VACF for 
the granular fluid described by (j2-ip and (j2-2p . Judging from the universality in the 
long-wave behaviors of granular flows and traffic flows, our result is expected to be 
robust and applicable to any situation of flows of dissipative particles. On the other 
hand, our result seems to be contradicted to the suggestion by Orpe and Kudrolli.'^ 
Although we expect that the interpretation of their result may be consistent with 
ours if the flow speed is small enough, we will need more systematic and careful 
treatments to clarify the relation between our analysis and their result. 

In this paper, we only discussed the linearly stable case of the fluctuations around 
the uniform flow (j3Tp . However, one of the most important characteristics of the 
optimal velocity fluid model (j2-ip and (j2-2p exists in the emergence of the congestion 
or separation between congested region and sparse region. In this sense, our analysis 
is far from completion. It is known that there is the frequency spectrum obeying 
y-4/3 £qj. granular flows ,13' 'I2Q),[2IJ ^j^g relation between the present analysis and the 
previous ones is not clear. We will have to clarify the relation in future. 

In conclusion, VACF in the granular flow governed by the optimal velocity fluid 
model ()2-ip and (j2-2p does not have the long-time tail but has the exponential decay. 
The amplitude of the exponential function becomes zero when there is no optimal 
speed Uq. Our result strongly suggests that we need careful interpretation of the 
experiment by Orpe and Kudrolli.'^ 

This work is partially supported by Ministry of Education, Culture, Science and 
Technology (MEXT), Japan (Grant No. 18540371). 

Appendix A 

The calculations of Cx{t) and Cz{t) 

In this appendix, we explicitly calculate ()3-9p and (j3-10p . For simplicity, we 
restrict our interest to 3-dimensional cases. 

For the calculation of Cx{t) and Cz{t) we need to know the expressions of cj+, 
(o"-(- — o"-) —1 and {((t+ — cr-)((T+ — cx)}^-*^ in the long- wave expansion as 

a+ ~ -z(co + C/o)t2A:^ + t ( eg - — ) kl - —kl + ■■■ (A-1) 
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2Tt 



m 



1 



(cj+ - o-_)(cr+ - o-_l) 



~ + Sicor'^A;^ + r'^ r 3 lOc; 



T 

I 

m 



2u k 



3Tt 



m 



2v\ki + 



(A-2) 



(A-3) 



Let us calculate Cx{t) at first. Substituting (jA-ip - (jA-3p into (j3-10p with some 
manupulation we obtain 



2^2 



cs) 



247r3m2 



dkk'^ exp 



-icokzt — ak t — { D -\ k<t 

m J 



(A-4) 



Introducing the cylidrical coordinate kx = k± cos 6 with f dk = f^^ dkj_ f^^ dO dz 
we obtain (j3-lip . where we use dA;2 exp[— (zcq/cz + aA;^)t] = ^J^^J(a£)e~'^^*/^'^'^'' 
and dkk^eM-(ikH] = {2/3t)-K 

Next, we calculate Cz{t). Substituting (jA-ip - (jA-3p into (|3-10p we obtain 



Cz{t) 



247r3 



m 



(A-5) 



where 



dkkz exp 



icokzt + akz + [ D + 



Tt 



m 



k'i \t 



■ exp 



2a3/2/3t3/2 
exp 



4a 



(A-6) 



4P + Zl)a5/2t3/2 



exp 



4a 



(A-7) 



where we use dzz'^ cos{cozt)e = ^?/2^3/2*"' c '^o*/(4a)_ Substituting (|A-6P and 
()A-7p into (jA-Sp we obtain (|3-12p as the long time behavior. 
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